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THERMODYNAMIC  THEORY  OF  ELASTIC  RELAXATION 
by  Josef  Meixner 

It 

The  Institute  of  Theoretical  Physics  of  Hheinisch-Westfalische 
Technische  Hochschule,  Aachen. 

(2.  Naturforschg.  9a,  654-663  [1954]}  received  on  May  11th,  1954). 


Assuming  that  the  relaxation  phenomena  of  a  system  may  be 
described  by  the  introduction  of  internal  variables  and  the  appli¬ 
cation  of  the  thermodynamics  of  irreversible  processes,  general 
properties  of  relaxing  systems  with  arbitrarily  many  external 
variables  are  examined.  In  particular,  the  transition  to  the 
after-effect  theory  of  relaxation  phenomena  is  carried  out,  several 
properties  of  the  after-effect  matrix  and  its  Laplace  transforms 
are  stated  and  proved,  the  thermodynamic  potentials  are  expressed 
in  terms  of  the  after-effect  functions  and  two  general  theorems  on 
relaxation  spectra  are  proved.  The  discussion  is  carried  through 
for  the  case  of  elastic  relaxation  taking  into  account  the  most 
general  properties  of  viscoelastic  bodies;  the  discussion  can 
immediately  be  extended  to  dielectric  relaxation. 

1.  General  Remarks  on  After-effect  Theory. 

The  after-effect  theory  of  relaxation  phenomena  reached  a  certain  stage 

1  2  12 
in  the  works  of  Gross  and  Hiedemann  and  Spence  (in  and  the  earlier 

literature  on  this  problem  is  surveyed)  for  the  most  important  case ,  namely 
that  of  the  isothermal  or  adiabatic  after-effect  between  two  conjugate 
thermodynamic  variables.  As  for  the  mechanisms  of  the  molecular  phenomena 
at  the  root  of  the  after-effect ,  there  are  few  quantitative  results  even 
in  cases  where  the  mechanisms  have  been  clarified  to  some  extent.  However , 
even  if  one  neglects  the  nature  of  the  mechanisms  and  examines  the  after¬ 
effect  theorv  only  in  its  phenomenological  form,  there  are  still  many  open 
questions.  Some  more  recent  work  does  deal  with  the  tensorial  formulation 


1  B.  Gross,  Mathematical  Structure  of  the  Theories  of  Viscoelasticity, 
Paris  ,  1953. 

2  E.  Hiedenar.n  and  R.  D.  Snence ,  Z.  Physik  133,  109[1952]. 
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of  the  after-effect  theory,  for  example  Sips'3  and,  in  particular,  Takizawa* 

5  6 

in  connection  with  the  research  of  Frenkel  and  Obrastzov  and  of  Oshida  ; 
however,  they  restrict  themselves  to  isotropic  bodies  and  only  Takizawa  drops 
the  adiabatic  condition  or  the  isothermal  condition;  he  also  includes  in  his 
consideration  the  temperature  after-effect,  thereby  taking  a  first  step  in 
the  direction  of  a  thermodynamical  after-effect  theory.  However,  what  is 
still  missing  is,  1)  the  formulation  of  a  thermodynamical  theory  which, in 
addition  to  the  after-effect  of  the  temperature  or  entropy,  includes  the 
after-effect  of  stress  and  strain;  2)  the  investigation  of  the  properties 
of  the  after-effect  matrix,  which  here  replaces  a  single  after-effect 
function;  3)  the  characterization  of  the  energy  conditions  within  the 
framework  of  the  after-effect  theory  (a  first  step  in  this  direction  has 
been  taken  by  Stavermann  and  Schwarsl7);  4)  the  replacement  of  the  after¬ 
effect  standpoint  by  the  introduction  of  suitable  internal  variables  and 
the  interpretation  of  the  after-effect  as  the  result  of  taking  limits  with 
the  internal  variables  tending  to  their  equilibrium  values;  and  finally 
5)  the  investigation  of  the  tensorial  properties  of  these  internal  variables, 

u 

in  particular  in  anisotropic  crystals.  Finkelstejn  and  Fastov  have  made  an 
important  contribution  to  the  last  problem  for  isotropic  bodies  by  the  intro¬ 
duction  of  the  relaxation  tensor. 

A  series  of  these  problems  is  dealt  with  in  this  paper  for  a  special 
case  of  the  general  after-effect  theory,  namely  for  the  case  which  is 
concerned  with  the  thermodynamics  of  irreversible  processes.  The  theory 
corresponding  to  this  is  the  so-called  thermodynamic  theory  of  relaxation 

3  R.  Sips,  J.  Polymer  Sci.  7,  19l[l951j. 

4  E.  I.  Takizawa,  Mem.  Fac.  Engng. ,  Nagoya  Univ.  5,  1[1953], 

5  I.  Frenkel  and  J.  Obrastzov,  J.  Rhys.  Chem.  USSR  2,  131[1940], 

6  I.  Oshida,  Mem.  Fac.  Engng.,  Nagoya  Univ.  2,  29[1950]. 

7  A.  J.  Stavermann  and  F.  Schwarzl,  Proc.  Roy.  Soc. ,  Amsterdam  B55, 

474,  486[ 1952  ]. 

8  B.  N.  Plnkelstejn  and  N.  S.  Fastov,  Dokl.  Akad.  Nauk.  SSSR, 
n.S.  71,  875[1950] . 
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phenomena  *  .  It  starts  from  the  assumptions  1)  that  the  thermodynamic 

state  of  the  material  concerned  may  be  represented  approximately  by  the 
internal  energy  as  a  function  of  the  entropy,  of  the  strain  components  and 
of  certain  internal  variables  which  describe  the  possible  internal  changes 
(chemical  reactions  in  the  most  general  sense)}  2)  that  the  usual  thermo¬ 
dynamic  formalism  yields  conjugate  variables ,  specifically  the  temperature  , 
the  stress  components  and  the  fluxes  corresponding  to  the  internal  variables 
and  finally  3)  that  the  variation  of  the  internal  variables  with  respect  to 
time  are  linear  homogeneous  functions  of  the  fluxes. 

All  the  considerations  can  be  carried  over  without  difficulty  to  the 
electrical  after-effect}  to  do  this  it  is  only  necessary  to  re- interpret 
the  symbols  and  formulae. 

2.  Thermodynamics  of  the  Anisotropic  Elastic  Body. 

We  consider  next  an  anisotropic  body  in  thermodynamical  equilibrium} 
this  can  be  unconstrained  or  constrained  or  also  in  the  sense  of  Schottky^ 
imaginarily  constrained.  We  denote  the  components  of  the  stress  tensor 

<rzz,  °2X>  axy>  ^  the  series  *1*  a2»  CT3> 
and  the  components  of  the  corresponding  strain  tensor 

£  %/£  >  £  %/S.  Hence  we  have  approximately  £  /V2  =6^  =  l/2(dv/dz  +  3w/dy) , 
if  v,  w  are  the  y-  and  z-components  of  the  displacement.  With  this  notation 
the  quantities 

6  6  6 

Z  CTi  *1'  Z^®*'  Z^  (1) 

i=  1  i=  1  i=l 

are  invariant  under  every  rotation  of  the  x,  y,  z  coordinate  system. 

The  state  of  the  body  is  determined  if  we  are  given  1)  the  Oj_  and 
the  temperature  T,  or  2)  the  C1  and  the  entropy  s  per  unit  volume,  or 


9  J.  Meixner,  Z.  Naturforschg.  4a,  594[1949). 

10  J.  Meixner,  Kolloid-Z  134,  3[1953]. 

11  W.  Schottky,  Thermodynamik ,  Berlin  1929,  p.  65. 
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3)  the  £  and  T  or  4)  the  a.  and  s.  As  reference  state,  we  denote  the 

*  4. 

stress-  and  strain-free  state  by  the  temperature  T  ;  it  can  also  be 

characterized  by  the  corresponding  entropy  s+  per  unit  volume.  We  always 

assume  that  all  the  variations  from  the  reference  state  i.e.  all  £  ,  a . , 

s  -  s  ,  T  -  T  are  so  small  that  linear  relations  exist  between  these 

quantities.  We  get  a  compact  notation  if  we  denote  s  -  s  by  £  and  T  -  T 

by  cr  .  Then  all  the  c\.  (i  =0,  1,  ... ,  6)  are  intensive  variables  and  all 
°  i  1 

the  £  (i  =  0,  1.,  . . . ,  6)  are  specific  (i.e.  referred  here  to  the  unit 
volume)  values  of  extensive  variables. 

For  the  internal  energy  u ,  the  free  energy  f  and  the  free  enthalpy  g , 
at  any  given  time  referred  to  the  unit  volume,  the  following  differential 


relations  are  valid 


du  =  Tds  + 


VJ 

£ 


de1  =  T+d£°  + 


£  -i 
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df  =  -sdT  +  oi  d  £x  =  -s+  dao  +  ^  d  E1  -  £°  d aQ,  (3) 

i=  1  i=l 


dg  =  -sdT  -  ^  £i  da^  =  - 


+  - 

s  do0  - 


e1  da±. 


In  what  follows,  unless  otherwise  indicated,  the  following  summation 
convention  is  now  observed:  -  any  latin  index  occurring  twice  in  a  product 
is  Wbe  summed  from  1  to  6  or  from  0  to  6  respectively  according  as  to 
whether  it  has  the  upper  index  '  or  not;  for  the  Greek  indices  which  arise 
later,  summation  is  similarly  over  all  internal  changes. 

In  the  domain  of  validity  of  Hooke's  law  extended  to  anisotropic  bodies 
and  with ,' moreover ,  small  variations  from  the  reference  state,  the  internal 
energy  can  be  expanded  to  a  first  approximation  as 

u  -  u+  =  t+ e°  +  1/2  (5) 


where  cik  =  cki  (i,  k  =  0,  1,  6)  (6) 

Hence  we  have,  from  (2),  the  linear  relations 

'Ji  =  elkEk  Ci-0.1 . 6)  (7) 
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ik  ki 

with  their  converse  (c  =  c  is  the  inverse  matrix  of  c^) 

t  —  c  °k  (i  =  0»  1#  ...  >  6).  (8) 

For  brevity,  we  call  c^  the  thermoelastic  matrix. 

There  then  results  for  the  free  energy  and  the  free  enthalpy 


f+  = 


+  ,  1 
-s  Oq  +  £ 


(ci’k* 


'oi,cok 


:'vi’  k»  J.ci'o  P i’ 
-It  £  +- - £  an 


oo 


00 


srz  °o°o> 
00 


(9) 


g 


+ 

g 


.ik 


CTi°rk* 


(10) 


4  -f  4* 

The  quantities  u  ,  f  ,  g  are  here  the  internal  energy,  free  energy  and 
f5ree  enthalpy  of  the  reference  state. 

In  (8)  we  have  assumed  that  the  61  are  uniquely  determined  from  the 
cr  and  conversely  (i.e.  that  the  determinants  of  the  matrices  c^  and  c 
are  different  from  zero  ).  Then  it  follows,  in  addition,  that  the  quadratic 
forms 

ik 

cik  xixk  and  c  xixk 

in  the  real  variables  are  positive  definite.  To  prove  this ,  we  start 
from  the  fact  that  the  internal  energy  of  the  adiabatically  strained  body 
with  £ 0  -  0  arid  the  free  energy  of  the  isothermally  strained  body  with 
<t  =  0  are  greater  than  the  internal  energy  and  the  free  energy  respectively 
of  the  reference  state.  Hence  the  quadratic  forms 

J  cok\ 

Cf,^,  xit  xk,  and  (c^.^i - -  -  — — ;  x^,  x^, 

in  six  variables  are  positive  definite  and  it  follows  that 


cik  xixk=  (ci'k' 


Coi,cok 


00 


-)  (°ok  xk^2>0’ 


if  the  x^  are  real  and  not  all  equal  to  zero;  for  is  the  inverse  thermal 
capacity  per  unit  volume  under  constant  strain  and  this  is  positive  because 
of  the  internal  thermodynamical  stability.  A  further  statement  of  the 
internal  thermodynamic al  stability  conditions,  which  imply  c^  x^x^  >0  for 


=  x. 


=  X, 


=  X, 


=  0, 


x0,  x1  =  real,  xQ  + 


#0 
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does  net  require  to  be  taken  into  consideration  separately*  since  it  only 
describes  a  special  case  of  the  positive  definite  character  of  the  quadratic 
form  which  has  already  been  proved.  In  the  same  way  as  c^  belongs 

to  a  positive  definite  quadratic  form  so  also  does  c^  by  known  theorems. 


3.  Consideration  of  Internal  Changes. 


We  now  extend  our  investigations  to  an  elastic  body*  in  which  internal 
changes  are  possible  and  whose  state  is  therefore  not  only  described  by  the 
G*-  or  by  the  cr^  or  a  certain  choice  of  seven  of  the  quantities  C1*  on  with 
different  indices  *  but  whose  internal  state  we  have  still  to  characterize 
by  so-called  internal  variables.  These  can  be  the  fluxes  of  the  internal 
changes  or  concentrations  of  independent  components  or  any  linearly  inde¬ 
pendent  combinations  of  these  quantities.  It  simplifies  the  notation  if 
we  assume  the  reference  state  *  now  also  to  be  characterized  with  respect 
to  the  internal  variables  *  to  be  stress-  and  strain- free  and  in  perfect 
thermodynamic  equilibrium  with  respect  to  the  internal  changes;  of  course* 
this  requires  it  to  be  perfect  only  with  respect  to  changes  which  can  take 
place  inside  the  period  of  observation;  other  changes  can  be  regarded  as 
not  permissible ,  their  fluxes  are  taken  into  the  thermodynamic  potentials 
and  relations  as  constant  with  respect  to  time.  In  considering  the  internal 
variables  we  shall  also  limit  ourselves  to  the  neighbourhood  of  equilibrium 
and  we  denote  by  the  deviations  of  the  internal  variables  from  the 
equilibrium  values  corresponding  to  the  reference  temperature  T+;  from  now 
on  we  call  these  differences  'f,a  the  internal  variables.  If  the  number  of 
internal  changes  is  n>  the  index  a  then  runs  from  1  to  n.  The  case  n  =  oo 
is  not  excluded;  it  is  then  assumed  that  the  infinite  series  which  occurs 
is  convergent. 

Then  for  the  internal  energy  it  follows,  generalizing  (5)  and  restricting 

terms  to  those  which  are  at  most  of  second  order  of  smallness*  with  c.,  =  c.  . 

ik  ki 

and  rap  =  rpa> 


u  - 


4*  4-  e  o 

u  =  T  G 


*1 


'ik 


Gi£k 


+  i  r 

2  o 


(id 


The  coefficients  c^*  q^a*  rflp  depend  only  on  the  reference  state,  i.e. 
T+  or  s+  respectively.  In  place  of  the  differential  relation  (2)  there 
holds  more  generally 


du  =  T+  d  G  0  +  ff  d  &  -  Aa  d£a 


on 

now 


(12) 
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where  the  coefficients  Affl  are  called  fluxes.  Then  from  (11)  and  (12)  follow 
the  linear  relations 


«i  -  “Ik  e“  -  "ia*‘ 


Aa  =  -V£  -rap? 


(i  =  0*  1  >  •  •  •  9  6)* 

^  (a  =  1,  2,  ,  n). 


(13) 


For  given  values  £*%  one  finds  the  equilibrium  values  of  the  'Z3a  by  minimissing 
the  internal  energy  or  equivalently  by  taking  the  fluxes  Aq  to  be  zero.  Since 
the  equilibrium  value  of  the  internal  energy  (i.e.  its  minimum  for  variation 
of  the  from  the  result  of  the  last  section  for  £  =  0)  is  greater  than  u  , 

except  if  all  £ 1  vanish ,  then  the  quadratic  form  in  the  7  +  n  variables  £T 


******  +  \ 


2  ik 


°P 


(14) 


is  first  of  all  for  £°  =  0  positive  definite;  from  analogous  properties  of 
the  free  energy  one  then  concludes ,  as  in  the  preceding  section ,  that  the 
quadratic  form  (14)  in  all  the  variables  G1  and  ^  is  positive  definite. 

cik  now  (^en0^es  the  thermoelastic  matrix  with  49  components  for  hidden 
coordinates,  i.e.  for  held  constant.  In  the  thermodynamic  equilibrium 
with  A  =  0,  it  follows  from  (11)  by  elimination  of  the  using  (13)  that 


u  =  T'  £u  +  i?(c.,  -  q. 

2  ik  Mia 


(15) 


if  r  p  is  the  reciprocal  matrix  of  r 


<xp* 


The  thermoelastic  matrix  for 


unconstrained  equilibrium  therefore  reads  c^  -  q^fl  r  "q^a  ;  also  it  belongs 
to  a  positive  definite  quadratic  form,  and  for  viscoelastic  bodies  (see 
below)  to  a  positive  semi-definite  quadratic  form. 

Now  there  are  some  particular  cases  still  to  be  investigated.  If 
det  r  „  =  0 ,  then  one  can  introduce  as  new  internal  variables  homogeneous 

I  CXP  I  n 

linear  functions  of  the  original  variables  '(£’ ,  in  such  a  way  that  rfflp  £ 
transforms  into  a  quadratic  form  with  fewer  than  n  variables;  also 
q.  £*  then  contains  no  more  of  the  same  new  internal  variables  than 

The  quadratic  form  (14)  is  then  only  semi-definite,  but  after 


va  -p 

rap>  >  4 

re-arranging  to  the  reduced  number  of  new  internal  variables  is  again  posi¬ 
tive  definite.  We  therefore  do  not  restrict  the  generality  if  we  assume 
(14)  positive  definite.  This  means  that  the  internal  variables  may  always 
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be  decomposed  by  a  suitable  choice  into  some  which  go  into  a  positive  definite 
quadratic  form  g  rQp  ^  in  u  ■  u  and  others  which  generally  do  not 


appear  in  u  -  u 

We  get  a  borderline  case  if  we  allow  the  body  being  considered  to  exhibit 

the  property  of  viscoelasticity.  Then  the  internal  energy  in  the  unconstrained 

thermodynamic  equilibrium  can  be  the  same  for  different  deformations j  thus  the 

internal  energy  of  the  isotropic ,  viscoelastic  body  in  equilibrium  does  not 

12  3 

depend  on  the  shear  deformation  ,  but  only  on  the  volume,  i.e.  on  £  +  £  +£  . 

The  quadratic  form  of  the  matrix  c^  -  is  then  semi-definite  and 

its  determinant  is  equal  to  zeroj  thus  this  case  can  exist  even  for  the  matrix 
Cik*  However,  in  the  following,  we  want  to  consider  exclusively  those  bodies 
which  in  constrained  equilibrium  (i.e.  for  constant  "E,0  or,  expressed  in 
another  way,  for  fast  stresses,  possess  no  viscous  behavior,  but  have  purely 


elastic  properties 

ik  ' 

of  the  matrix  c  ,  the  reciprocal  of  c^,  with  det 

Besides  the  internal  energy  (11)  we  still  estimate 


£  0  and  the  existence 


Hq 


Then  we  can  always  assume  det  c . 

T4  0. 

l 

the  free  energy  and 
the  free  enthalpy,  as  they  can  be  calculated  from  their  definitions  using 
(13),  These  are,  in  suitable  variables, 


f  - f+  ■  -s+%  - 


2c“  °o°o 
00 


ci'o 


o  o 
00 


^ci’k ' 


coo 


c— 

00 


(q. 


(16) 

00  r  00 


q~,q. 


ia 


g  -  g+  =  -s\  -  |  ^  a  A  *  I  O’) 

“1th  o/  -  rflp  =  rap  -  q^c11^  .  (18) 

Finally,  £°  can  be  calculated  from  (11)  as  a  function  of  u  -  u+, 

lc  *  Q 

£  to  an  equivalent  approximation,  i.e.  including  the  quadratic  terms 

e>°  1  (■>,  +>  C°°  t  +\2  Cok*  /  +uk'  ci'k*  _i'  vi 

£  =  -  (u  -  u  ) - xr  (u  -  u  )  -  -j;.-  (u  -  u  X  -  - x-  e  £ 

T+  2T+3  T+2  T+ 


^oa  /  +\_a  ^i'a  <4’, a  raB  _a_B 

-5*?(u-u  >5  -  —  ^5  V- 


(19) 
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Using  the  differential  relations  (2),  (3)  and  (4),  one  again  obtains  the 
linear  relations  (13)  between  both  sets  of  conjugate  variables  in  a  different 
form.  Then  it  follows  from  (17)  that 


C1  =  crk  -  qQ1  |a(i  =  0,1,  ...  ,  6)  , 

(20) 

Aa  =  q£CTk  -  ?a^P(a  =  lj  2»  n)j 

+ 

while  (19)  leads  exactly  to  the  equations  with  —  ~  —  in  place  of  £? . 

T 


Now  the  thermodynamic  relations  can  be  completed  using  the  relations 
which  describe  the  change  with  respect  to  time  of  the  internal  variables 
%a .  Since  the  d  T^a/dt  vanish  for  A^  =  0 ,  i.e. ,  for  thermodynamic 
equilibrium,  then  we  can  obviously  express  the  a/dtj  at  least  in  the 
neighbourhood  of  equilibrium  (which  is  always  assumed  here) ,  as  linear 
functions  of  the  A 

a 

=  E  ^Ap  (a  =  1,  2,  •••)  n).  (21) 


This  substitution  gives  exactly  the  thermodynamics  of  irreversible  processes' 
and  it  says  moreover  that  the  Onsager’s  reciprocity  relations 

E®*3  =  EP“  (22) 


are  valid. 

Also  the  matrix  Ea^  belongs  to  a  positive  definite  quadratic  form,  if 


the  local  entropy  generation;  however,  if  it  were  semi-definite,  then  its 

determinant  would  vanish  and  it  would  give  a  linear  combination  of  the  2-a , 

which  remains  constant  with  respect  to  time  and  can  be  different  from  zero. 

Hence  we  can  also  write  (21),  by  introducing  the  matrix  E  „  which  is  the 
aB 

reciprocal  of  E  p  and  with  Eap  =  Epa  •  as  follows 


(a 


1,  2,  . . .  ,  n). 


(23) 


*12 


12  J.  Meixner,  Z.  phys.  Chemie  B  53,  253[1943];  Ann  Phys.  (5)  43, 
244[  1943  ]. 
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4.  Transition  to  the  After-effect  Interpretation. 

The  change  with  respect  to  time  of  the  state  of  the  material  considered 
depends  in  the  thermodynamic  theory  of  relaxation  phenomena,  which  has  just 
been  discussed,  only  on  the  instantaneous  state,  i.e.  on  the  totality  of 
independent  variables  describing  the  state.  Therefore  these  are  here  perhaps 
the  C1  and  the^°.  The  change  with  respect  to  time  of  the  conjugate  quanti- 

w 

ties  <7^  is  certainly  no  longer  determined  if  one  does  not  know  the  %  ,  unless 
one  knows  other  measures  of  the  system  besides  the  o  at  the  instant  concerned 
the  at  earlier  instants  would  be  suitable.  We  get  just  such  a  represen¬ 
tation  of  the  if  we  eliminate  the  in  the  relations  of  the  last  section. 
For  this  purpose  we  start  from  the  2nd  equation  in  (13)  and  eliminate  now 
the  ka  using  (23).  We  obtain,  with  the  notation  d/dt  =  D, 

(r»p  *  W>)5fl  =‘<Ikae1'-  <24> 

We  denote  by  (D )  the  operator  matrix  which  is  the  reciprocal  of 
rQp  +  by  substituting  for  from  (24)  in  the  first  equation  (13), 

there  results 

ai  =  ^Cik  ~  qia#  ^D^\p  ^  £  *  (25) 

The  structure  of  this  equation  becomes  considerably  clearer  if  we  make  use 

of  the  fact  that  both  the  symmetric  matrices  r  „  and  E_D  with  positive 

ap  ap  | 

definite  quadratic  form  can  be  transformed  by  a  real  matrix  T^“  into 
diagonal  form.  Therefore  let 

T;Sr +  v,Te- =  (1 +,t«  D):6«e 


[1  *  t0D 


for  a  =  £ 


for  a  ^  £ 


(26) 


be  universally  valid  in  D  with  positive 
for  constant  £^.  Then  it  follows  that 


*ap(D) 


the  so-called  relaxation  times 


1  +  TyD 


(27) 
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and  with  the  abbreviation 


from  (25) 


’ia 


cr. 

1 


QiV  QkY  1 

i  +  t  d 

r 


(28) 


(29) 


This  transformation  to  principal  axes  therefore  effects  the  partial  fraction 
decomposition  of  the  operator  #a^(D).  The  relation  (29)  becomes ,  after 
multiplication  by  the  denominators  1  +  t^D,  a  dynamic  thermoelastic  relation 

between  the  a.  and  the  o  together  with  their  derivatives  up  to  nth  order. 

1  13 

It  is  an  extension  of  the  well-known  dynamical  equation  of  state  of  acoustics  . 

From  (29)  we  now  come  immediately  to  the  relationship  between  the  and 

the  in  the  sense  of  the  after-effect  theory.  For  any  function  x(t) ,  it  is 

true  that 


r  co 


x(t)  =  (1  +tD) 


-v/t 


x(t  -  v)dv 


(30) 


and  her.  e  ye  can  transform  ^29)  to 


n  00 


T«rf  Trll 


*  cik e  (t)  • 


(v>V  Q.v  0 

ik  /  lDr  kt 

Y 


cik(v)  £k(t  -  v)dv. 


-V/Ty 


(31) 

(32) 


First  of  all  note  that  (29)  is  a  system  of  inhomogeneous  differential 
equations  for  the  cn  (t) ,  as  one  recognizes  by  multiplying  out  the  denominators 
1  +  n^D.  (3D.)  :is  only  a  particular  integral  of  this  inhomogeneous  system  of 
equations?  however,  it  describes  the  behavior  of  the  system  of  interest  to 
us,  because  one  may  always  regard  the  solutions  which  are  decreasing  exponentials 
with  respect  to  time  as  having  died  out,  if  one  has  left  the  system  to  itself 
for  long  enough  before  the  start  of  an  experiment. 


13  L.  I.  Mandelstam  and  M.  A.  Leontovitsch ,  J.  Exp.  Theor.  Phys.  U.S.S.R 
7,  438[D  '"37]. 
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Meantime  we  start  from  equation  (20)  instead  of  from  equation  (13)  and 
again  eliminate  the  5°,  and  obtain  analogously 


where 


n  oo 


e1(t)  = 


C^CT.  (t)+ 


C^v) 


ok(t 


v)dv, 


QiYQkY  £ 


-v/cY 

~ 


(33) 


(34) 


j_v  Y 

with  real  quantities  Q  and  positive  relaxation  times  t'  for  constant. 
Therefore  we  have  here  as  the  antithesis  of  the  stress  relaxation  for  constant 
strain  [expressed  by  the  minus  sign  in  (32)]  a  retarded  strain  for  constant 
stress.  For  this  reason  one  frequently  calls  the  ty  retardation  times. 


5.  The  Properties  of  the  After-effect  Matrix. 

We  consider  in  addition  to  the  after-effect  matrices  C^(v)  and  C^k(v) , 
their  Laplace  transforms 


r  oo 


-pv 


*  r  ^iY  ^Y 

Cik(v)dv  =  — - - 1  , 

1  +tyP 


(35) 


= rki(p)  = 


roo 


e-pv  Cik(v)dv  =  ££ 
0  1  +  xYp 


(36) 


If  one  starts  with  =  0,  t1  =  0  for  t<0  and  constructs  the  Laplace 
transforms  of  equations  (31)  and  (33),  one  recognizes  that  both  the  matrices 
c..  -  P-.  (p)  and  c*k  +  P^(p)  are  reciprocals  of  one  another.  Hence  we 


Jik 

have 


[cik  “  rik(p)][ckl  +  rKX(p)  ]  = 


kl, 


(37) 


The  parametric  value  p  =  0  corresponds  to  thermodynamic  equilibrium ,  and 
the  parametric  value  p  =  oo  to  the  constrained  equilibrium  -  realized 
during  very  fast  changes  of  the  c  (or  of  the  <y^  respectively),  in  particular 
for  periodic  stress  of  very  high  frequency.  The  thermoelastic  matrix  for 
thermodynamic  equilibrium  is  therefore  c^  -  Pjjf(0);  it  has  a  non-negative 
quadratic  form.  In  particular,  it  is  positive  definite  if  in  the  equilibrium 
not  only  the  cr.  are  uniquely  determined  by  the  G1 ,  but  also  the  inverse  is 


-  13  - 


true.  It  is,  however,  positive  semi-definite  if  the  C1  are  no  longer 
uniquely  given  by  the  then  det  c±k  -  P^O)  |  =  0.  Thus  this  relation 
characterizes  in  general  the  viscoelastic  behavior  and  indeed  does  so  not 
only  for  isotropic  materials,  but  also  for  single  crystals.  For  isotropic 
materials,  however,  it  follows  in  particular  from  the  vanishing  of  this 
determinant  that  all  third  order  minors  also  vanish  and  hence ,  that  an 
isotropic  material  with  this  property  for  any  shear  deformation  is  visco¬ 
elastic. 

The  functions  H^Cp)  and  P1^(p)  are  regular  analytic  functions 
except  at  p  =  -l/r  and  p  =  -1/ty  respectively.  For  the  matrices  P..  (p) 
and  P  (p)  one  easily  obtains  from  (35)  and  (36)  the  following  statements: 

For  real  xi  and  'fa,  p>  0  the  quadratic  forms  THP^p)*^  and  '7&'Pik(p)xixk 
are  non-negative  definite.  Under  the  same  conditions  the  quadratic  forms 
%,[c  +  Tik(p)]  xixjc>  ”  r^(p)  Jxixk  are  positive  definite,  however 

for  the  last  form  with  the  restriction  p  ^  0,  if  the  body  is  viscoelastic. 
Finally,  it  is  true  for  real  x^,  which  are  not  all  zero  and  )&p.>0  that 

[Cik  -rilc(p)]xixk>7^|  [o^  -  r^O)]  >0;^,p[c:Uc  +Pik(p)>0.  (38) 

We  shall  return  to  the  meaning  of  the  last  two  statements ,  which  are  also 
true  for  materials  whose  after-effect  properties  no  longer  come  under  the 
thermodynamic  theory  of  relaxation,  in  another  connection. 

The  relationship  between  the  £ 1  and  the  cu  is  closely  analogous  to  the 
relationship  between  the  in- flowing  currents  in  an  electric  2m-pole  (m  =  7) 
and  the  applied  voltages.  If  we  associate  the  with  the  voltages  and  the 
d£ /dt  with  the  intensities  of  the  currents  of  the  2m-pole  then  the  impedance 
matrix  of  the  electrical  network  corresponds  exactly  to  the  matrix 
^ik(P)  =  ^  ^Cik  Therefore,  if  it  is  possible  to  realize  the 

matrix  Z^^(p)  by  a  suitable  network,  we  obtain  an  electrical  model  for  the 
whole  elastic  and  thermal  behavior  of  our  body.  For  the  case  m  =  1  (with 
t’-'o  conjugate  variables)  electrical  models  with  2-poles,  in  particular, 

RC-models  with  resistances  and  capacitances ,  have  long  been  used  as  a  means 
of  visualizing  relaxation  rhenomena. 

We  observe  further  that  the  after-effect  functions  (^(v)  and  CT'(v) 
for  real  v>0  belong  to  non-negative  definite  quadratic  forms;  this  follows 
immediately  from  (32)  and  (34).  Even  more  strictly,  it  is  true  that  these 
quadratic  forms 
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Cik^xiXk 


and 


Cik(v)xixk, 


(39) 


for  those  real  x^  for  which  they  are  not  identically  zero  in  v,  are  strictly 
inonotonic ally  decreasing  for  0<  r<  ooj  i.e.  their  nth  differential 
coefficient  is  different  from  zero  and  has  algebraic  sign  (-l)n.  This 
property  is  characteristic  of  the  structure  of  the  after-effect  functions. 
From  it *  by  a  theorem  on  strictly  monotonic  functions ,  it  follows  that  the 
quadratic  form  C^v^x^.  may  be  described  as  a  Stielt jes  integral 


poo 


-W 


d[<ik(x)xixk 


h 


(40) 


with  the  function  ^jc(^)xix]c  n°t  monotonically  decreasing  in\.  A  similar 
result  is  true  for  C^(v).  Let  us  restrict  ourselves  -  something  that  is 
not  too  serious  and  corresponds  to  the  case  of  discrete  relaxation  spectrum  - 
to  functions  ^^(X)  >  which  apart  from  ju:mp  discontinuities  =  l/z^>0 
are  constant*  then*  at  each  jump  discontinuity* 


^ik^Y  +  "  °^XiXk 


is  non-negative  definite*  therefore 


“'ik 


R. 


ik,Y 


(41) 


where  R^.^x^x^  is  a  non-negative  definite  quadratic  form  for  every  Y. 

Then,  however,  R^  ^  may  be  represented  as  a  sum  of  at  most  s  +  1  (<  7) 
products 

Rik,Y  =  RiY  K?  +  RiY  ^Y  +  *•*  +  RiY  Kf  (i'  k  =  O’1’*  •  *  »6)  (42) 


with  real  factors  *  and  if  one  counts  the  relaxation  time  t  as  s  +  1  -  fold  * 
one  returns  directly  to  (32). 

6.  The  Thermodynamic  Potentials  in  the  After-effect  Interpretation. 

First  of  all  we  eliminate  the  '£a  in  (11)  and  obtain 

U  -  U*  =  T*  e°  4  \  <>  eV  .  e1  ek  +  \  (-^il —  ,  <«) 

*  ^  1+tD  *  1  +  t  D  1+tD 

Y  Y  Y 
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whence  ,  by  (37)  and  (32),  it  follows  further  that 


u(t)  -  u*  »  T+e°(t)  +  g  cikei(t)^(t)  -  £i(t)  Cik(v)£k(t  -  v)dv 


oo  r  oo 


0  JO 


C*  (v  +  w)£i(t  -  v)ek(t  -  w)dvdw  ; 

XK 


where  C^v  +  w)  means  the  derivative  with  respect  to  the  argument. 
Correspondingly  we  obtain  from  (17) 

g(t)  -  g+  =  -s+cro(t)  -  |  oi(t)ok(t)  +  a±{  t)  Cik(v)ck(t  -  v)dv 

J  o 


oo  r»co 


0  JO 


(v  +  w)o-i(t  -  v)ak(t  -  w)dvdw 


Farther  we  have 


poo 

T+[s(t)  -  s+]  =  u(t)  -  u+  -  J  gi(t)  £k(t)  +  e^t)  C±k(v )  (t-v)dv 

JO 


pco  poo 

’Jo  l 


C^k(v  +  w)  £^(t  -  v)  £k(t  -  w)dvdw, 


if  we  put 


5°  =  B  “JIL  ,  g1  =  8i(i  =1,2,  ...  ,  6).  (47) 

T 

An  analogous  representation  may  be  obtained  for  the  free  energy. 

Fran  these  results  it  therefore  follows  that  the  thermodynamic  potentials 
at  time  t  may  also  be  represented  by  the  values  of  their  independent  variables 
at  time  t  and  at  earlier  times  and  also  that  only  the  after-effect  matrices 
enter  into  them.  The  same  statement  is  true  for  the  local  entropy  generation 


roonco 

0  0  «*<»♦.) 

JU  JU 


dv  dw 


'T00  ik  d<j.(t  -  v)  dak(t  -  w) 


dv  dw* 
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7.  Theorems  on  Relaxation  Spectra. 

If,  after  some  preceding  strain  and  entropy  changes,  we  hold  all  the 
£*  constant  with  respect  to  time  from  t  =  0  onwards ,  then  from  t  =  0  onwards 
the  £a  given  by  (24)  adjust  themselves  to  equilibrium  values.  Then  for 
t  >  0  every  can  be  written  as  a  sum  of  exponential  functions  with  exponents 
-t/x  ,  apart  from  a  constant  (the  term  containing  S1).  We  call  the  t  the 

u  t  j  G 

relaxation  times  for  constant  Z  and  their  totality  we  call  the  relaxation 
spectrum  Sp(£°,  ...,  £6)  for  constant  £^j  we  count  each  relaxation  time  n 
times,  where  n  =  its  multiplicity.  In  the  dynamic  equation  of  state  (29) 
not  all  the  need  occur j  for,  all  Q^(i  =  ••• »  6)  can  vanish  for 

certain  Y.  We  want  to  denote  the  t^.,  actually  appearing  in  (29)  and 
different  from  one  another,  as  the  effective  relaxation  spectrum,  i.e. 
operative  for  the  thermcelastic  relations  (29),  Spe££.(£°,  ...,£®)  for 
constant  £.  If  we  ask  further  how  a  definite  one  of  the  cr . ,  perhaps  aQ> 
approaches  the  equilibrium  value  for  constant  £  ,  then  under  the  circum¬ 
stances  again  only  a  subset  of  the  effective  relaxation  times  is  of 

significance;  we  call  them  briefly  the  q-effective  relaxation  times  and 

lc 

their  totality  the  q-effective  relaxation  spectrum  for  constant  £  ,  and  use 
the  notation  Sp0^(qj  £°,£^,  ...,'-^)  for  it. 

One  can  now  define  relaxation  spectra  for  other  conditions,  while  one 
holds  constant,  instead  of  all  ZT  >  only  some  of  them  and,  in  addition,  the 
fqr  the  remaining  indices.  Such  a  choice  we  call  a  'selection*  of 
variables.  As'  an  example,  we  therefore  investigate  the  spectrum  for  constant 
£°»  £i,  ... ,  £q,  ,  ... ,  o>  and  call  it  the  spectrum  for  this  selection  of 

variables.  Hence  we  solve  equation  (13)  for  cr^cr^,  •••  >  ,  ... ,  P° 

and  obtain 

q  6 

t7i=Z1ik£k*  2  +  o,i, 

k=  0  k=  q+1 

q  6 

£i  =  -  X  +  Z  Vic  +  ♦la  =  q  +  X>  "•  >  6)’ 

k  =  0  k=q+l 

q  6 

Ao  =  -  ^  4a£k  +  X  *kaCTk  “  Pap  ~  1»  2>  *  n^*  (49) 

k=  0  k=q+l 
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The  summation  convention  is  now  valid  only  with  respect  to  a.  The  relations 

^"ik  =  ^ki*  nik  =  nki  and  p0p  =  pf!a 3  as  wel1  33  *he  chosen  signs  follow 
immediately  from  the  fact  that  the  liniar  relations  (49)  can  be  derived  from 
a  potential  \Jr  with  the  differential 

q  6 

= y  ^  af-k  -  Z eV  v*? 

k=  0  k=q+l 

Rirther ,  it  follows  immediately  from  the  evaluation  of  u  =  u  +  ~  cr.C1- 

1  &  Ct 

that  lilc»  n^  and  pa^  belong  to  positive  definite  quadratic  forms.  In 
addition,  we  sill  have  the  relations 


V  $- ■  A«(«  - a* 


•  •  •  > 


n). 


(50) 


By  a  suitable  linear  transformation  of  the  ,  we  now  restore  to  the 

form  6  Q  and  E  Q  to  the  form  t  6  Q.  The  V  are  then  obviously  the  relaxation 
op  ap  §  op  a  J 

times  for  constants  ~ ,  ...,£,  cr  •••»  If  we  assume  on  the  other 
hand  £°,  ...,  cr^,  ...,  cr6  constant,  then  by  (49)  and  (50) 


qq 

-with  no  loss  in  generality,  we  put  £  0  =  ...  =  =  a  =  ...  c  =  o,  which 

is  equivalent  to  replacing  the  by  the  differences  of  the  "v  and 
their  equilibrium  values  for  constants  C°,  ...,  erg*  The  solutions  of  this 
system  of  equations  are  again  sums  of  exponential  functions  with  exponents 
~t/ra  and  one  shows,  just  as  in  ®  taking  into  account  1  >0,  that  by 

ordering  the'  xa  and  the  xa  hy  magnitude 


T1^T1*£  T2  ^n^n 


(51) 


is  true.  We  write  briefly  for  this 
Sp(£  ,  •  •  • ,  £^  ,  £^>  ^q+i*  •  •  •  5  cjg )^  ®p(£  >  , 


aq*  aq+l* 


.  •  • , 


ag> 


and,  in  particular,  write  here  <  instead  of^if  in  (51)  the  equality  sign 
is  nowhere  valid.  In  words,  we  say  that  the  first  spectrum  is  'not  larger' 
than  the  second,  or  in  the  other  case  that  the  first  spectrum  is  'smaller* 
than  the  second.  Note  that  the  <:  -  or  orelation  defined  in  this  way  is 

Spj  are  true ,  then  the 


not  transitive}  in  fact,  if  Sp^<  Sp^  and  Sp^ 
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relaxation  times  of  the  first  spectrum  are  certainly  smaller  than  those  of 

the  third  spectrum,  however  they  need  not  be  interlaced.  Since  the  special 

choice  of  variables  to  be  held  constant  is  obviously  immaterial  for  this 

result,  we  can  state  the  following  in  general 

If  one  replaces  a  specific  extensive  variable  £q  by  the  corresponding 

intensive  variable  cl ,  in  a  selection  of  variables  e1,?,  (the  totality  of 

q .  k 

the  indicies  of  the  £  and  c*k  therefore  contains  each  of  the  numbers  from  0 

to  6  exactly  once) ,  then  the  relaxation  spectrum  for  the  system  is  not 

diminished. 

We  now  consider  further  the  two  spectra 


and 


Speff(QjE  ° 9  el>  eq,<Vl» 

«  /  _,o  *  1  ^q-1  \ 

Speff(qj  £,£,...,£  ,  oq>  ...»  a6) 


By  eliminating  the  £>a  from  (49)  and  (50)  with  £° 


we  obtain 


i  fS  \ 

*qq  'qq 

1  +  ^„D 
°  / 


.=  <r6=  0, 
(52) 


The  right  side  can  be  written  as  a  rational  function  in  D  with  as  the 
degree  of  denominator  and  numerator,  if  is  the  number  of  the  xa  actually 
appearing  in  (52)  and  different  from  each  other.  While  the  xa  make  up  the 
spectrum 

Speff(qj  e°»  e1*  •••>  eq?  <vi»  ...»a6) 

one  obtains  the  spectrum 

Speff(q*  e0,el>  eQ  1'  °V  a6)t 


if  one  looks  for  the  zeros  of  the  numerator  of  the  rational  function  mentioned 

and  puts  them  equal  to  the  negative  reciprocal  relaxation  times  of  the  latter 

spectrum.  FTom  the  behavior  of  the  rational  function  for  -l/ra  and  because 

1  -it  0,  it  follows  that 
qq  Pqprqa 

3peff (qj  £°,.«.  »  £q»  cfq^i >•  •  • »  £  >**«>£q  »  <Tq,#**  9 

The  effective  relaxation  times  of  these  two  spectra  are  therefore  different 
from  each  other  in  that  between  every  two  times  of  the  one  spectrum  lies 
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one  of  the  other  and  the  smallest  relaxation  time  of  the  first  spectrum  is 
smaller  than  that  of  the  other.  Also  we  can  generalize  this  result  as 
follows 

If  one  replaces  a  specific  extensive  variable  E^  by  a  corresponding 

intensive  variable  o_  in  a  selection  of  variables  £■*",  <r.  ,  then  the  q-effective 

q  k 

relaxation  spectrum  is  enlarged. 

If,  in  particular,  1  -  =  0  is  true,  then  a  relaxation  time  of 

the  second  spectrum  is  infinitely  large  and  we  have  viscoelasticity. 

The  two  theorems  obtained  are  independent  of,  and  thereby  valid  for, 
whichever  spatial  coordinate  system  the  E1  are  referred  to.  However, 
in  the  proof' we  have  never  employed  the  fact  that  the  o\  ,  E1  for  i  ^  0  are 
components  of  the  stress-  and  strain-tensors ,  but  only  made  use  of  the  fact 
that  the  linear  relations  (13)  may  be  derived  from  a  positive  definite 
quadratic  form  (14)  over  the  differential  relations  (12).  Therefore,  the 
two  theorems  on  relaxation  spectra  are  also  valid  if,  instead  of  the  cr ,  £  , 

_  if  1 

we  use  any  linear  combinations  cr.  ,  -  u!  cr,  and,  as  conjugate  quantities, 

k  i  k  *k^*  I  i  I 

the  c  determined  by  E1  =  u^.  t  ,  where  det  u^1  ^  0  has  been  assumed. 

This  fact  is  important,  for  example,  if  one  considers  the  relaxation  between 

pressure  and  volume.  Then  it  depends  on  the  connection  between  cr,  +  or0  +  <r- 
12  3  x  a  o 

and  £  +  £  +  E  perhaps  for  constant  shear  deformations  or  constant  shear 

stresses  with,  in  addition,  constant  temperature  or  entropy.  It  follows  on 

putting 

f1  =  —  (E1  +  £2  +  £3) ,  £2  =  -i-  (E1  +  £2  -  2£?) ,  £3  =  —  (E1  -  Z2) 

vr  y/r  vr 


and 


-  2c r,),  a,  = 


y/r 


a2) 


and  retaining  the  other  cr^,  &  with  the  indices  i  =  0,  4,  5,  6.  Then,  as 
•selections’  of  variables  £° ,  ,  €2 ,  . ..,  £6  and  £°,  a^,  . ..,  £6  are 

to  be  chosen  for  adiabatic  conditions.  For  isothermal  conditions ,  in  both 
selections  £°  is  to  be  replaced  by  aQ> 

These  statements  on  the  relaxation  spectra  are  valid  in  the  same  way 
for  discrete  relaxation  spectra  with  infinitely  many  relaxation  times ,  and 
are  also  valid  for  continuous  relaxation  spectra.  Prom  the  relation  (51) 
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it  follows  in  particular  that  ,  if  the  spectrum  of  the  t0  fills  a  finite 
continuous  interval  a<  xa<;  b ,  the  spectrum  of  the  xa  fills  the  same 
continuous  interval}  there  is  however ,  at  least  in  those  cases  where  one 
deals  with  a  q-effective  spectrum,  as  one  sees  in  particular  examples 
(see  '*'0)  ,  still  an  isolated  relaxation  time  a  outside  this  interval. 

A  corresponding  statement  is  true,  if  the  x  fills  a  continuous  interval 
a<  Ta<:  b;  then  there  is  an  isolated  relaxation  time  T>b. 

The  statements  of  this  section  are  closely  connected  with  the  principle 
of  Le  Chatelier-Braun. 

8.  Periodic  Changes  of  State 

If  the  Z1  are  periodic  functions  of  the  time  with  time  dependence 
e^t ,  then  it  follows  from  (29)  that 


ai  = 


V.V.  V, 

1  +  ilOT  J 

r 


(53) 


i.e.  the  connection  between  the  cr^  and  C1  is  formally  the  same  as  in  the 
equilibrium  («  =  0) ,  but  with  a  frequency  dependent  complex  thermoelastic 
matrix. 

For  slow  changes  of  state,  for  which 


1  for  s  =  m  +  1,  m  +  2,  ... , 
s 

one  can  write  as  an  approximation 


n 


m 


i  |  [file  2  QiYQkr]  "  X 


ScY 


Y=  m+1 


Y=  1  1  +  3WTY 


00 


2 


Y=  m  +  1 


V  V  Tt 


l  6* 


Among  these  assumptions  the  fastest  relaxation  phenomena  have  been  regarded 

as  having  practically  died  out!  they  make  themselves  felt  only  in  a  formal 

n 

change  of  the  thermoelastic  matrix  c^->-  c^-  V1  -  it  corresponds 

Y=  m+1 

to  a  thermodynamic  equilibrium  for  which  the  fast  relaxation  processes  have 
died  away  and  the  slow  are  continuing  -  and  in  an  extra  term  of  viscous 
character  proportional  to  it.  Here  it  is  not  a  question  of  a  proper  viscosity, 

V, 

as  one  sees  if  one  solves  (53)  for  the  £  and  then  again  approximates  for  the 
fast  relaxation  processes  in  the  same  way  as  above;  in  fact  there  is  no  term 


of  the  sum  with  the  factor  On  the  contrary,  proper  viscous  behavior 

(viscoelasticity)  exists,  as  we  have  already  noted,  only  if 
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